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ABSTRACT 


A FORTRAN computer code for the reduction and analysis of experimental heat transfer data has been 
developed. This code can be utilized to determine heat transfer rates from surface temperature measurements made 
using either thin-film resistance gages or coaxial surface thermocouples. Both an analytical and a numerical finite- 
volume heat transfer model are implemented in this code . The analytical solution is based on a one-dimensional, 
semi-infinite wall thickness model with the approximation of constant substrate thermal properties , which is 
empirically corrected for the effects of variable thermal properties . The finite-volume solution is based on a one- 
dimensional, implicit discretization . The finite-volume model directly incorporates the effects of variable substrate 
thermal properties and does not require the semi-infinite wall thickness approximation used in the analytical model. 
This model also includes the option of a multiple-layer substrate. Fast, accurate results can be obtained using either 
method . This code has been used to reduce several sets of aerodynamic heating data, of which samples are included in 
this report. 


INTRODUCTION 


Information on aerodynamic heating rates is a critical factor in the design of hypersonic vehicles. A 
common experimental technique for the determination of aerodynamic heat transfer rates is wind tunnel testing of 
models which have been instrumented with surface temperature sensors such as thin-film resistance gages or coaxial 
surface thermocouples. These sensors are used to measure surface temperatures during a wind-tunnel test, and the 
surface heat transfer rates are then calculated from the recorded time histories of the surface temperatures. 

An analytical relationship between the heat transfer rate and the measured surface temperature-time history 
was developed by Vidal (ref. 1), and Schultz and Jones (ref. 2). This relationship is based on the model of one- 
dimensional heat transfer to a semi-infinite substrate which has constant thermal properties. Cook (ref. 3) extended 
this theory to include substrate materials with variable thermal conductivity. Hartunian and Varwig (ref. 4) and 
Miller (ref. 5) developed approximate correction factors to account for substrate materials with variable thermal 
conductivity and thermal diffusivity. The one-dimensional heat transfer problem with variable thermal properties has 
also been solved numerically by many authors, for example White (ref. 6), Patankar (ref. 7) and Dunn et al (ref. 8). 
Pittman and Brinkley (ref. 9) and Bradley and Throckmorton (ref. 10) have extended the numerical method to 
multiple-layer substrate problems. 

In this report, the One-Dimensional Hypersonic Experimental Aero-Thermodynamic data reduction code, 
1DHEAT, will be detailed. 1DHEAT is an fast, user-friendly FORTRAN code for heat transfer data reduction which 
incorporates both analytical and numerical models. Variable thermal properties are accounted for in both models, and 
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both models can be used to reduce either thin-film gage data or coaxial thermocouple data. Additionally, multiple- 
layer substrates can be dealt with using the finite-volume model. 


SYMBOLS 


Ch Stanton number (dimensionless heat transfer coefficient), q / [pu„ (h^ - )] 

c p specific heat (J/kg-°K) 

E sensor voltage (mV or V) 

h enthalpy (J/kg) 

H dimensional heat transfer coefficient q / (h^ -h^) 

k thermal conductivity (W/m-°K) 

L( ) Laplace transform of ( ) 

q heat transfer rate (W/m 2 ) 

Q heat addition (J/m 2 ) 

R sensor resistance (£1) 

r adiabatic wall recovery factor 

RMS( ) root-mean square of ( ) 

S( ) standard deviation of ( ) 

s Laplace transform variable 

T temperature (*K) 

t time (sec) 

u velocity (m/s) 

x depth through substrate(m) 

x* dimensionless penetration depth 


a 

«, 

p 

p' 

<i> 

X 


thermal diffusivity, k / ( pc p ) (m 2 /s) 
coefficient of resistance ( 1/°R) 
thermal product, Jpc^k (W-s 1/2 /m 2 -”K) 
correction factor 

? k 

transformed temperature, 0=1 — dO (°K) 

o k ° 

dummy time integration variable (s) 
mesh Fourier number, <r = 

(Ax) 2 


2 



e 


temperature defined by 9(x y t) = T(x,t) - T(x, 0) (*K) 
density (kg/m 3 ) 


P 


SUBSCRIPTS 


0 

1 

2 

amb 

cal 

gage 

int 

I 

r 

IP 

ref 

s 

tot 


freestream 

initial (pre-test) condition 

start of averaging window 

end of averaging window 

ambient (pre-test) 

calibration condition 

surface temperature measured by gage 

cell interface 

left boundary 

right boundary 

reference point temperature 

reference condition 

surface temperature 

stagnation 


ONE-DIMENSIONAL HEAT TRANSFER 


Figure 1 depicts a generic surface temperature sensor mounted on a substrate material. For a thin-film gage, 
the sensor is typically a platinum or palladium film, and the substrate is a thermally insulative material such as 
Macor or quartz. For a coaxial thermocouple, the sensing element is a thin junction of two thermocouple materials, 
such as chromel and constantan, while the substrate is the body of the thermocouple. In either case, the sensor is 
typically designed so that the thickness of the sensing element is much less than that of the substrate. The sensing 
element, therefore, has a negligible effect on the heat transfer to the substrate, and the temperature measured by the 
sensing element is identical to the temperature at the surface of the substrate. It is also assumed that there is no 
lateral conduction through the substrate and that heat is conducted only in the direction normal to the surface. The 
final assumptions are that the substrate is of infinite depth, and that the temperature rise at infinity is zero. 
Together, these assumptions are referred to as the one-dimensional, semi-infinite solid model. Using this model, the 
temperature distribution in the substrate can be given by the partial differential equation 
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dx\dx) p dt 


(la) 


The thermal properties (thermal conductivity, specific heat and density) are all functions of temperature. 
The boundary conditions are given by 


T(x,t = t 0 ) = T 0 

T(x = 0j) = T s (t) (lb) 

q(x = ~ 0 = 0 


A common simplification is the assumption of constant thermal properties, in which case the temperature 
distribution is given by 


d 2 T \ dT 

dx 2 a 0 dt 

where 

k 

a 0 = = constant 

P c p 

The boundary conditions are the same as for equation (1). 


TEST-TIME AND THE SEMI-INFINITE ASSUMPTION 

The semi-infinite substrate assumption, which is required for an analytical (but not numerical) solution of 
equation (1), is only approximately correct during an actual experiment. The accuracy of the assumption is a 
function of material depth, diffusivity and test duration, and can be measured in terms of the increase in the 
temperature or heat transfer rate at the back of the substrate. The relevant dimensionless parameter for the thermal 
penetration depth is given in reference 2 


* 2 Vat 


Expressions for the back face heat transfer and temperature for a constant heating rate were also derived in 
reference 2. 


Ll 

T. 


J-* r 


(*' '4n}erfc(x ’) 


(4a) 
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q. 


(4b) 


= erfc(x') 


The period of the semi-infinite test time can be estimated from Figure 2 (which is redrawn from reference 2) 
as a function of the thickness of the substrate. Ideally, the substrate should be thick enough that the ratio qjq, at 

the back face of the substrate be no more than a few percent for the test time desired. Because the penetration depth 
is inversely proportional to the square root of diffusivity, materials with lower diffusivity (such as ceramics) have a 
longer semi-infinite test time, whereas high diffusivity materials (such as metals) have a shorter test time. 


HEAT TRANSFER MODELS 


Equation (1) can be solved to determine the temperature distribution, and thus the heat flux, either 
analytically or numerically (through the use of the finite-volume technique). Analytical solutions are possible for 
materials in which the thermal properties are constant or in which only the thermal conductivity varies with 
temperature; an empirical variable properties correction to the constant properties solution is also possible. Variable 
thermal properties are treated directly in the finite-volume solution method. These methods are summarized in this 
section, and are discussed in detail in Appendices A and B 

Constant thermal properties 

The simplest analytical solution (ref. 2) is found by assuming that the material properties are all 
independent of temperature. The surface heat flux can then be found from the measured time history of the surface 
temperature by 




T: ~ T, 


(5) 


where the thermal product 


Po = 4P c p k ( 6 ) 

is evaluated at the ambient temperature. The derivation of equation (5) is presented in Appendix A. 

Alternatively, the total heat added to the substrate can first be calculated using the method given by Kendall 
and Dixon (ref. 11) 
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■At 


Qn = 


_/*o 


Tt + Ti. i 




(7a) 


and the heat transfer rate can then be computed with the finite-difference approximation used by Hedlund et al (ref. 

12 ) 


q( t ^-q _ dQ„ 2Q n _ 8 - g n _ 4 + g n+4 + 2Q n+g 


dt 


40 A t 


(7b) 


The schemes represented by equations (5) and (7) are referred to as the direct method and the indirect method, 
respectively. In general, the instantaneous heat transfer rates computed using the two methods are not identical. 
This is because the temperature difference term in the numerator of equation (5) of the direct method tends to 
accentuate fluctuations in q, whereas the temperature sum term in the numerator of equation (7a) and the wide 
differencing stencil in equation (7b) of the indirect method tend to smooth fluctuations in q . However, over a given 
time interval the average values of the heat flux will be approximately equal. The indirect method then, would be 
preferable in a situation where the random noise level is relatively high in comparison to the heat transfer rate. If 
however, the fluctuations in q were of interest, as for example in a turbulent flow, then the direct method would be 

preferable. Further discussion and comparisons between these two methods are presented in the Data Reduction and 
Analysis section. 


Variable thermal conductivity 

Equations (5) and (7) are based on the assumption of constant substrate thermal properties. In hypersonic 
heat transfer experiments, the increase in the temperature of the substrate is usually not negligible. This assumption 
then is not valid, and the above formulations are valuable only as starting points for more detailed analyses. 

A more accurate analytical solution can be found by assuming that the thermal conductivity of the substrate 
is a function of temperature but that the thermal diffusivity is a constant (ref. 3). 

k 

a = constant 

PCp ( 8 ) 

This is a reasonable approximation for metallic substrates, such as that of a coaxial thermocouple. With 
this assumption, the equations for the direct and indirect methods become 

Direct method; 


M t ) - V Z ft-i 


(9) 
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Indirect method: 


where 


°-Skl 


4>, + <P,-i 




-A/ 


(10a) 


q«») = q n =^ 

at 


- 2 Q,-8-Qn-4+Q„ + 4+2<2^ s 

40 At 


(10b) 
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Variable thermal conductivity and thermal diffusivity 

Equations (9) and (10) are inadequate for the insulative substrates typically used with thin-film gages, such 
as Macor, quartz, or Pyrex. For these substrates, both the thermal conductivity and thermal diffusivity are functions 
of temperature. For this case, there are no exact closed form solutions for heat transfer such as given by equations 
(5,7) or (9,10). However, an accurate approximation for the heat transfer rate may be made through the use of a 
correction factor for the variation of material properties with temperature. This approach was first used by Hartunian 
and Varwig (ref. 4) for quartz and Pyrex substrates. Correction factors were also derived for Pyrex by Cook (refs. 
13,14) and for Macor by Miller (ref. 5). New correction factors have been derived for Macor and quartz by 

comparison with finite-volume solutions, as detailed in Appendix C. The formulation for heat transfer with variable 
material properties is 


where 


fy(T) - ip, H + P'AT S ] ( 12a ) 

AT s = T g«ge ~ T amb (12b) 


where /?' is the correction factor which accounts for the variation of material properties and q fio is the constant 
material properties heat transfer rate computed using either the direct (eq. 5) or indirect method (eq. 7). 

Finite-volume 

When all thermal properties are functions of temperature, an "exact" solution to equation (1) can only be 
obtained through a finite-volume model of the substrate. 
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T i = T(x„n 


(13) 


The numerical model is more flexible than the analytical model in that the semi-infinite assumption does 
not need to be made; this assumption is replaced with a physical boundary condition at the interior (back-face) of the 
substrate. The boundary condition may be that of an adiabatic (zero heat transfer) back face, a constant temperature 
back-face, or a temperature measurement made using a second sensor located at the back face of the substrate. This 
finite wall thickness analysis is most useful because it is not uncommon for a wind tunnel model to have a wall 
thickness small enough that the semi-infinite assumption is violated during a test, which would result in incorrect 
heating rates being computed from a semi-infinite analysis. The finite-volume analysis also makes it possible to 
model a multiple-layer substrate. The finite-volume method is detailed in Appendix B. 


1DHEAT CODE DESCRIPTION 


OVERVIEW 

The analytical and computational models discussed in the previous section have been incorporated into 
1DHEAT. 1DHEAT is a flexible, user-friendly FORTRAN code for the reduction and analysis of experimental heat 
transfer data from thin-film and coaxial thermocouple surface temperature sensors. 1DHEAT has been compiled and 
executed on SUN and Hewlett-Packard workstations. Execution time is dependent on the number of temperature 
sensors, length of the data acquisition window, and the heat transfer data reduction model employed, but is typically 
on the order of a minute. As an example, a test with data from 70 thin-film gages with 4 seconds of data taken at a 

50 Hz sampling rate (14000 data points) has been reduced in under 30 seconds on both Hewlett-Packard and SUN 
workstations 


INSTALLATION AND SETUP 

The 1DHEAT code and sample data sets are distributed in a tar file , '1DFILES', which is accompanied by 
an installation script file, TDINST'. To install the 1DHEAT code, copy the files TDFILES' and 'lDINST to the 
host computer and run the installation utility by typing "lDINST". The TDINST' utility creates a home directory 
1DHEAT' for 1DHEAT and installs the source code there. It also creates subdirectories for the sample data sets and 
installs the necessary raw data files and input files therein. The source code 'lDHEAT.f may then be compiled and 
immediately be used to reduce and analyze the sample data sets, as will be discussed in the section "Data Reduction 
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and Analysis . To reduce new data sets, raw data input files and setup files must first be created, as will be detailed 
in the next sections. 


NAMING CONVENTION AND DIRECTORY TREE 

The 1DHEAT code uses or creates several files for each data set which is analyzed. In order to avoid 
confusion, 1DHEAT data files' names follow a specific naming convention, and are organized in an hierarchical 
directory tree based on wind tunnel test series and run number identifiers. The 1DHEAT directory tree is depicted in 
Figure 3. The naming convention for 1DHEAT data files is TD lestIDrunID .ejrf'. where test ID ' and ninlD ' are 
the identifiers for the test series and run number identifier of a data set. These identifiers may each be any 
combination of up to six numerals and/or characters (upper or lowercase). The extension 'ext' denotes the file types, 
which are summarized in Table 1. A file name example is TDT293R005.ex/', which would designate the 5th run in 
test series 293. 

For each test series, the user must create a subdirectory 'tes ttestID ' under the parent directory in which the 
1DHEAT source code resides. For each run in a test series, the user must create a subdirectory 'run runID ' in the 
test testID ' subdirectory. All input files for a given run must be placed in that run subdirectory; 1DHEAT 
automatically places the output files in the appropriate ran subdirectories. 


FILE DESCRIPTIONS 

Input files 
TDtestlDrunlD.inp' 

TD testIDrunID fvinp' 

' 1 D testIDrunID .de gk' 

TD testlDrunTD .vnlt' 

To reduce experimental data with the 1DHEAT code, a setup file and a raw data file for each data set which 
is to be analyzed must be provided by the user; if the finite- volume model is to be used in the data reduction, an 
additional finite-volume setup file is also required. 

The setup file, TDtestlDrunlD.inp’, contains entries on flow properties, gage and substrate types and 
other information which is required for the reduction of the raw data. The finite-volume setup file, 
1 D t estlDrunip . fvinp' contains information on the number of substrate layers, the thickness of each layer, and the 
number of computational nodes to be used for each layer. Examples of each of these file types are given in 
Appendix D along with descriptions of the entries in the files. 
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The raw input data for 1DHEAT may be either the sensor output voltage-time histories or the temperature- 
time histories computed by the user from the voltage-temperature calibrations of the sensors. The raw data files for 
these options are entitled TDtgsHPrunIP.volt’ and TD testlDrunlP .degk'. respectively. If raw voltages are to 
be input, the voltage-temperature calibrations must be of the same form as those in the 1DHEAT code. For thin- 
film gages, the voltage-temperature relationship used is based on the calibration methods used at the Langley 
Research Center: 


T = 


amb 


AE 

^ amb a R 


-tl + a K (T amb -T„ l y] 


(14) 


where subscript cal' refers to calibration conditions, subscript 'amb' refers to ambient pre-test conditions and A £ is 
the measured gage output (relative to pre-test baseline) in millivolts. For coaxial gages the conversion is the NBS 
standard for Type-E (chromel-constantan) thermocouples 

T = T r P + ^a i E‘ (15) 

1=0 

where E is the absolute measured gage output (which must be in volts for the NBS curve fit) and 7^ is the reference 
point temperature of 273 *K (freezing point of water). The coefficients of equation (21) are listed in Table 2. If 
gages with different voltage-temperature calibrations than these are utilized, the user must first convert the voltage- 
time history to a temperature-time history. 

Output files 
'lDtestIPrunID.t' 

'lD testIPrunIP .qx' 

TD testlPrunlP .chx' 

’ 1 DtestTPrunlP.distx' 

TD testlPrunlP. distxnew' 

' 1 D testlPrunlP. qresx' 

’ 1 D testIPrunIR .chresv' 

' 1 DtestIPrunIP . tback' 

' 1 DtestlPrunlP.g/ip^/Vf 

Several output files are created for the results computed using each of the data reduction methods. The 
method is identified by the number V in the file extension, which corresponds to computations based on: 
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x=l) Indirect Analytical Method: 

Equation (7) with the correction in (12) for thin film gages, or equation (10) for coaxial thermocouples. 



x=2) Direct Analytical Method: 

Equation (5) with the correction in (12) for thin film gages, or equation (9) for coaxial thermocouples. 
x=3) Finite-Volume Method: 

Numerical solution of equation (13) for either gage type. 

The output files , lD fest^DrunIP .qx , and T 1 D tes tID run ID .chx’ contain the time histories for each gage 
of the heat transfer rate and the Stanton number, respectively. The TD testIP^unIP .distx , file contains the 
averaged surface distributions of the heat transfer rate, q y Stanton number, C A , the dimensional heat transfer 
coefficient, //, and the extrapolated reference temperature heating rate, q refy as well as the standard deviations of the 
averaged heat transfer rate and Stanton number and their root-mean square values. The surface heating distributions 
are all based on the time-averaging window specified in the input file. The user may recompute these values based 
on a different averaging window, in which case the new distributions are written to the TD testlDrunlD .distxnew’ 
file. The output files TD lestIPrunID .qresx , and , l DtestIP^unIP .ch^esx , contain the time histories of the heat 
transfer rate residual and Stanton number residual, respectively. Finally, the TD testIDruTlID .t , file contains the 
temperature-time history for each gage. The quantities above are all defined in the Heat Transfer Computations 
section. 

The TD testIPrunIP .tback’ and TD test^P^unID .gggg/^ , files are unique to the finite-volume method. 
The TD testlPrunlD .tback 1 file contains the temperature-time history of the back face node of each gage, and can 
be used to estimate the validity of the semi-infinite assumption. The T DtestIPriinID .gfiffg/ff file contains the 
temperature profile through the substrate of sensor "gageid" at a user-specified time. 

Error log 

The 1DHEAT code generates an error log entided TDHEAT.errorlog’ which resides in the 1DHEAT root 
directory. This file contains a record of any errors which occurred during a session and is automatically output to the 
screen at the end of the session. The error log can be useful in identifying problems with the setup or data files, and 
as a record of any errors which occur if 1DHEAT is run in batch mode. 


UNITS 


All input data used by 1DHEAT must be in SI units except for the thin-film gage coefficient of resistance, 
a Ry (in the TDtestIPrunIP.inp , file) which has units of (l/’R) (to agree with the current format for gage 

calibration data at NASA LaRC) and the gage surface position data (also in the TD testIDrunIP .inp 1 file) which 
may be in any system of units. No computations are made with the position data; it is only included to allow the 
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user to plot surface heating distributions. All output data is in SI units with the exception of the gage position data 
in the T DlestIDrunID distx' and ' 1 DtestIPru nID .distxnew' files. 


HEAT TRANSFER COMPUTATIONS 


The 1DHEAT code computes the time histories of the heat transfer rate for each gage using the method 
specified by the user. After the heat transfer rates have been computed, the dimensional heat transfer coefficient, H, 
and the dimensionless heat transfer coefficient (Stanton number), Cy,, can be computed by 


H = 


£ 

(Ah^-AhJ 


and 


Q = 


(Ah^, — A/t*, ) 


(16) 


(17) 


where the enthalpies are defined by: 


M w = h(T w )-h(T rp ) (18a) 

M mv = h(T m ,)-h(T rp ) (18b) 


The factor hiT^) is the enthalpy at a user-defined reference point temperature. The enthalpies in equations 
(18a) and (18b) are written as differences to remind the user that the computations for the wall and adiabatic wall 
enthalpies MUST be based on the same reference point temperature. For perfect gas flows, the reference point 
enthalpy is usually taken as zero at 0' K, while for chemically reacting flows it is sometimes taken to be zero at the 
standard temperature of 298 < ’K. The reference point in 1DHEAT is 298‘K. 

The adiabatic wall enthalpy is related to the total enthalpy by 

Ah aw = Aha +(r -l)-y (19) 

and the recovery factor, r, is 1 for a stagnation point and for a flat plate is typically taken as Pr 12 for laminar flow 
or Pr 1/3 for turbulent flow. Note that in equation (19) it is assumed that the edge velocity is equal to the freestream 
velocity, and thus the Stanton number computations in 1DHEAT are based on freestream as opposed to boundary 
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layer edge conditions. This approximation is made in order to simplify the setup file formats and reduce the amount 
of data that the user must provide. 

A "reference" heating rate can also be computed. This is an estimate of the value of the heat flux at a 
specified wall temperature (such as the ambient temperature), and should not be confused with the reference 
temperature, T of boundary layer computations. Assuming that the heat transfer coefficient for a particular gage 
remains constant with time (allowing for experimental noise and measurement error), the reference heat transfer rate 
can be extrapolated from an averaged Stanton number by: 

‘W = C h ' (A /l w - (20) 

The recovery factor, reference temperature, freestream velocity and total enthalpy are all user inputs in the 
setup file, while the wall enthalpies are computed in 1DHEAT from the temperature-time histories of the sensors. 


STATISTICS 


After the heat transfer computations are made, 1DHEAT performs a statistical analysis of the heat transfer 
data by computing the RMS (root mean squares) values and the standard deviations of the averaged heat transfer rate 
and Stanton number. The RMS values are given by 


RMS(q) = 


hr- I fc) ! 


av 8 n=t, 


( 21 ) 


and the standard deviations are given by 


S(g) = 


"1 




n avg ^ 


( 22 ) 


The formulas for the Stanton number values have the same form. In equations (21-22), i,and (i, + n avg ) are the time 
indices of the start and end of the averaging window. 

The standard deviation is indicative of the experimental uncertainty of the time-averaged values, which can 
be influenced by factors such as flow quality and gage electrical noise. If the data is free from electrical noise 
produced either by the sensor or the data acquisition system, and the flow quality is high (and laminar) and the heat 
transfer rate (or Stanton number) is constant, then the standard deviation should tend toward zero. Data from 
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individual sensors which show significantly higher standard deviations than other nearby sensors should be considered 
suspect as this suggests poor sensor performance. Overall higher standard deviations may suggest either system-wide 
electrical noise or significant freestream unsteadiness or turbulence. 

The normalized time history of the residual of the heating rates has been found to a useful parameter with 
which to characterize the wake flow-establishment process in short-duration (test times of milliseconds or less) test 
facilities (Hollis and Perkins, ref. 15). The normalized residual is given by 


= 



<7; 


(23a) 


or for the Stanton number by 





(23b) 


The derivatives in equation (23) are numerically approximated by a four-point central derivative of the same 
form as equation (10b) for the indirect method, and by two-point central derivatives for the direct and finite-volume 
methods. The establishment process over the entire body is characterized by an RMS of the heat transfer rates of all 
the sensors (with a similar form for the Stanton number): 


RMS{hq) i = 


— — (Aip + A^-K. 
ngages v 





(24) 


Because the residual values are dependent on the heat transfer model and the numerical derivative 
approximation, they should be considered a qualitative, not quantitative source of data Also note that this RMS 
computation is spatial (based on data from gages at different positions at a given time) as opposed to the temporal 
computation (based on data from the same gage at different times) of equation (21). 

DATA REDUCTION AND ANALYSIS 

The data reduction and analysis options available in 1DHEAT are listed in the 1DHEAT main menu, which 
is shown in Figure 4. The options are divided into two groups: semi-infinite options and finite-volume options. 

Options 1 and 2 correspond to the semi-infinite solution methods of equations (5) and (7) with the 
correction for variable thermal properties of equation (12) for thin-film gages, and to equations (9) and (10) for 
coaxial thermocouples. Options 5 and 6 correspond to the finite-volume numerical method of equation (13). The 
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output from any of the above options consists of files containing the temperature-time histories 
( 1 D feSflP rpnip .t'), the heating rate-time histories (T DtestlDrunlD .qx'V and the Stanton number-time histories 
of each gage ('lDtestlDrunfP.chx'). Additionally, the averages and statistics for the heating computations over 
the specified time window are output to (TD testIDrunID .distx'V Options 3 or 7 are used to recompute these 
averages over a new time-averaging window (output to M ptestlDrunlD .distxnewT Options 4 or 8 are used to 
output the residual time histories of the heat rate and Stanton number. Finally, option 9 is used to compute the 

temperature distribution (using the finite-volume method) through the substrate of a single gage at a specified time 
(output to 'lD testlDrunlD .gggg/fiO- 

These options can be explored using the sample data sets created by the installation utility. The sample 
data is identified by the test series "demo"; in this test series are two runs: "coax" and "thin" (the file name templates 
for these data are thus lDdemocoax.ex/ and lDdemothin.gr/, respectively). The "democoax" data set is from a run in 
the NASA Langley Research Center 31" Mach 10 Air Tunnel (this facility is described by Micol in ref. 16) of a 2" 
radius hemisphere instrumented with Type-E coaxial surface thermocouples. The "demothin" data set is from a run 
in the same facility of a 70 sphere-cone (0.5 nose radius, 0.05" comer radius, and 2" base diam eter) with thin-film 
gages on a Macor substrate. 

A 1DHEAT session is illustrated on the following pages using the "democoax" data as an example. 
1DHEAT output is shown in bold face and user keyboard input in shown in italics. 

First, the raw data from the "democoax" set is reduced using the semi-infinite methods (direct and indirect): 

I " 1DHEAT" v2 .20 PROGRAM OPTIONS: | 

' i 

I SEMI- INFINITE METHOD | 

I 1 . Reduce raw data (volts) . | 

I 2 . Reduce raw data (degrees K) . | 

I 3 . Recompute statistics over new I 

I time interval. I 

I 4. Compute residuals. I 

! I 

I FINITE -VOLUME METHOD | 

I 5 . Reduce raw data (volts) . I 

I 6 . Reduce raw data (degrees K) . | 

i 7 . Recompute statistics over new | 

I time interval. I 

I 8. Compute residuals. I 

I 9 . Temperature profile for a gage at a | 

I specified time. i 

| I 

I 10. Quit . I 


Enter your choice. (1-10) 2 

Enter test ID (up to six characters) demo 
Enter run ID (up to six characters) coax 

Opened setup file; t estdemo /runcoax/lDdemocoax . inp 

Opened temperature data file: t estdemo /runcoax/ lDdemocoax . degk 

Reading data for 7 gages 
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Read 501 data points per gage for 7 gages 
Computing dQ/dt ( 351 points per gage),.. 

Integration window: 1.72 (sec) to 8.32 (sec) 

Finished computing dQ/dt 

Computing q ( 351 points per gage)... 

Integration window: 1.72 (sec) to 8.32 (sec) 

Finished computing q 

Computing statistics for 7 gages 

Averaging window: 4.0 0 (sec) to 5.00 (sec) 

Writing data... 

Finished writing data for 
Test: demo 

Run: coax 

After examination of the reduced data to determine a more exact data-averaging window, 1DHEAT is used to 
re-average the data over a new time interval: 


"lDHEAT" v2 . 2 0 PROGRAM OPTIONS: 

SEMI - INFINITE METHOD 

1 . Reduce raw data (volts) . 

2 . Reduce raw data (degrees K) . 

3 . Recompute statistics over new 
time interval. 

4. Compute residuals. 

FINITE -VOLUME METHOD 

5 . Reduce raw data (volts) . 

6 . Reduce raw data (degrees K) . 

7 . Recompute statistics over new 
time interval. 

8. Compute residuals. 

9 . Temperature profile for a gage at a 
specified time. 

10. Quit . 


Enter your choice. (1-10) 3 

Enter test ID (up to six characters) demo 

Enter run ID (up to six characters) coax 

Opened setup file: t estdemo /runcoax / lDdemocoax . inp 

Opened reduced data file: testdemo/runcoax/lDdemocoax . ql 

Opened reduced data file: testdemo/runcoax/lDdemocoax . q2 

Opened reduced data file: testdemo/runcoax/lDdemocoax . chi 

Opened reduced data file: testdemo/runcoax/lDdemocoax . ch2 

Reduced data runs from 1.72 (sec) 

to 8.32 (sec) 

Enter start time (sec) of averaging window. 3.7 5 
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Enter end time (sec) of averaging window. 4.25 
Computing statistics for 7 gages 

Averaging window: 3.75 (sec) to 4.25 (sec) 

Writing data 

Finished writing data for 
Test: demo 

Run: coax 

Additional data sets can now be reduced, or as in this example, the program can be exited. 


M 1DHEAT " v2 . 20 PROGRAM OPTIONS: 

SEMI -INFINITE METHOD 

1 . Reduce raw data (volts) . 

2 . Reduce raw data (degrees K) . 

3 . Recompute statistics over new 
time interval. 

4. Compute residuals. 

FINITE -VOLUME METHOD 

5 . Reduce raw data (volts) . 

S . Reduce raw data (degrees K) . 

7 . Recompute statistics over new 
time interval. 

8. Compute residuals. 

9 . Temperature profile for a gage at a 
specified time. 

1 0 . Quit . 


Enter your choice. (1-10) 10 

SESSION ENDS 


The reduced data from this session and from reduction of the "demothin" data set are useful in illustrating 
the features of the 1DHEAT code. Data from the "democoax" set is plotted in Figures 5 through 9. Data from the 
"demothin" set is plotted in Figures 10 through 17. Flow conditions for these two runs are given in Table 3. 

The measured temperature-time histories of the thermocouples on the hemisphere (seven, distributed along 
a single ray at 15’ increments) from the "democoax" data set are plotted in Figure 5. This data is found in the 
lDdemocoax.t file. Note that because of the high diffusivity and conductivity of the thermocouples, the 
temperature rise, and thus the increase in wall enthalpy, was small during the test. Because the wall enthalpy 
increase was small the heat transfer rate (which is proportional to the difference between the flow stagnation enthalpy 
and the wall enthalpy) remained nearly constant through the duration of the run. 

The time histories of the heat flux at the stagnation point gage (6 = 0) as computed using both the direct 
and the indirect method are shown in Figure 6. These time histories are found in the TDdemocoax.ql' and 
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’lDdemocoax.q2' files. Although the two time histories follow each other closely throughout the run, it can be seen 
that the direct method highlights fluctuations in the heating whereas the indirect method tends to smooth out these 
fluctuations. 

The stagnation point heating-time history is shown with that of the Stanton number (from the 
'lDdemocoax.ql' and TDdemocoax.chl’ files) in Figure 7. The Stanton number remains constant as is generally 
expected, and because of the small temperature rise during the run, so does the heat transfer rate. Note that the data 
shown in Figure 7 (and Figure 6) extends well beyond the semi-infinite test time estimated from Figure 2, which 
was approximately 1.5 seconds from the point at which heating began. Because the semi-infinite assumption was 
eventually violated, and because lateral conduction also began to become significant as the surface temperature 
increased, the one-dimensional analysis of the data became less accurate with time. This is why the Stanton number 
begins to increase slightly after the 4.5 second mark. In most cases (including this one), a sufficient amount of data 
can be acquired before the one-dimensional model is invalidated, and the later data are usually ignored. This portion 
of the run is shown here only to illustrate the limits of the applicability of the one-dimensional model. 

The average values for the "democoax" data set of the indirect and direct method heat fluxes (from the 
'lDdemocoax.distl' and TDdemocoax.dist2' files) measured over the hemisphere are plotted in Figure 8, and as 
would be expected from the time histories in Figure 6, are nearly identical. Finally, in Figure 9 the averaged direct 
method heat transfer rate values are shown along with their standard deviations (both in the TDdemocoax.dist2' file); 
the small values of the standard deviations indicate that the data from this run are quite reliable. 

The measured temperature-time histories of the thin-film gages on the 70° sphere-cone (twelve, distributed 
along a single ray at increments of S/R b = 0. 1) from the "demothin" data set are plotted in Figure 10. These data are 
found in the 'lDdemothin.t' file. Because the thin-film gages were mounted on an thermally insulalive substrate 
(Macor) they experienced a much greater temperature rise than the coaxial thermocouples. 

The time histories of the heat flux at the stagnation point gage (S/R b = 0) as computed using both the 
direct and the indirect method are shown in Figure 11. These time-histories are found in the 'lDdemothin.ql' and 
'lDdemothin.q2' files. Because of the large surface temperature rise the direct and indirect method heating-time 
histories do not remain constant as was the case for the coaxial thermocouples, but do again follow each other 
closely throughout the run. Because of this, an averaged heat transfer rate is a meaningless quantity. However, as 
shown in Figure 12, the Stanton number (from the 'lDdemothin.chr file) is nearly constant and can be averaged. 
To illustrate this point, consider the distributions plotted in Figures 13 and 14. Averaged heat transfer rates are 
plotted in Figure 13 and as expected, the average values changes with the window position. The average heat flux in 
the original window of 3.0 to 4.0 seconds (from 'lDdemothin.distl' file) is different than that in a user-input (using 
option 3) window from 2.75 to 3.25 seconds (from the 'lDdemothin.distl new') file. Both are different than the value 
extrapolated to a reference temperature of 300 °K (also in ’lDdemothin.distl’). On the other hand, the averaged 
Stanton number values do not change with the averaging window, as shown in Figure 14. From these figures, it 
can be seen that while both the Stanton number and heat flux computations at a given time are valid, it is usually 
simpler to refer to the Stanton number since the heat transfer rate only makes sense if information on the wall 
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temperature is provided (or if the wall temperature rise is small in comparison to the stagnation temperature, as seen 
in the coaxial thermocouple data set). 

The finite-volume data reduction options were also employed in the reduction of the 'demothin' data set. As 
shown in Figure 15, the direct and indirect method Stanton number-time histories, which are constant property 
solutions multiplied by an empirical correction factor, are in excellent agreement with the numerical solution (from 
the TDdemothin.q3' file) until well into the run when the semi-infinite approximation begins to become less 
accurate. 

Temperature profiles through the substrate of the stagnation-point gage ("TF19") from the finite-volume 
solution at various times during the run are shown in Figure 16 (these data are from the TDdemothin.TF19' file, 
which was created using option 9). In Figure 16, note that at 3.5 seconds, the temperature at the back-face has 
begun to rise. The back-face temperature is one indicator of the validity of the semi-infinite substrate model. When 
the back-face temperature rise begins to become significant (relative to the sensor temperature) the semi-infinite 
model does not apply, since a constant back-face temperature is assumed in this model. The back-face temperature- 
time histories for each gage are located in the TDdemothin.tback' file; these data are plotted for several gages in 
Figure 17. 


FINITE-VOLUME APPLICATIONS 

Use of the finite-volume options in the 1DHEAT code makes it possible to analyze problems for which the 
classical semi-infinite substrate model (indirect or direct method) is not applicable. Specifically, the 1DHEAT code 
can be used to reduce data from gages with substrates of finite thickness, and these substrates may be comprised of 
one of more layers of different materials. Although the limiting assumption of one-dimensional conduction still 
applies in these cases, the finite-volume approach represents the actual physics of the problem more accurately than 
the semi-infinite substrate model. 

In Figure 18, the heating distribution on a 70‘ sphere-cone similar to the one in the "demothin" data set is 
plotted for a test in the Langley 20" Mach 6 Air Tunnel (ref. 16). The distributions shown in Figure 18 are from 
analysis by the indirect, semi-infinite substrate method and by the finite-volume method. In this example (courtesy 
Thomas Horvath, NASA Langley), the model was 6” in diameter and was machined from chromel and instrumented 
with Type-E coaxial thermocouples. This model was designed for testing in a shock tube with run times on the 
order of milliseconds, as opposed to the ran times in the Mach 6 tunnel which were on the order of seconds. Testing 
in the Mach 6 tunnel led to the semi-infinite assumption being violated in the comer region of the model, where the 
chromel substrate was considerably thinner than for the rest of the model. Because of this, analysis by the semi- 
infinite method led to invalid results in the comer region. Results from the semi-infinite analysis showed that there 
was a broad local heating peak at the comer of the model, which was much greater than that of the 70* sphere-cone 
model in the demothin" example. However, when the finite-volume model, which accounts for non-infinite 
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sute thicfmess, was employed for daa raducdon, this peak was foam, h. considerably ^ which made 
distnbution more consistent with the "demothin" data set. 

““'"Pl-'ayer substrate results are illustrated to Figure 19 to which heating dara from two 0.25' radius 4" 

NAsItuT 7 'TT Uly “ ' he N “ A ^ 2 °" MaCl ’ 6 ^ ^ iS (-tmesy Scot, Berry 

g ey). e irst cylinder was machined from Macor and insttumented with thin-film gages A 0 002" 

: 0fUP ' J 7 ^ « *o machtoed from Mmor. This UpL lay 2s 

m frumeuted with thto-fUm gages (note the tighter gage spacing possible on the Upi.es, see the discus ton o^ 

Utoto to Appendix C). When tested, bod, cylinder were ahgned pnrpenditodar to the freesfream flow so toa, , 

.ton-film gages, which weie to a single row along the leading edge of each cylinder, measured the two-dimensiona, 
stagnation point heating along the length of each cylinder. 

thickness* 2 Uol'T” IT" “ ‘ hem ' a ' ^ ^ "" “ WaS " Ulch *»- «- toe 

P x ayer, which made this a mulnple-layer problem. Thus, the heat transfer was computed 
using the finite-volume technique, with 5 points to the Upilex layer (dx . 0.0,02 mm, and with 244 points in the 

rr ' ( " 7 60 “* 19 - uumputed St»o„ iTra oI I 

Upilex/Macor cylinder with that from toe all-Macor cylinder, fra which dato was .educed using the indirect semi- 

Lolls' “ me ' h0<i ' The C0, " P “ ,e ' "-ber, distohutions for Ihe two 

modeb which were geometncally identical but had dissimilar thermal response characteristics due lo the difference in 

U S e materials, are in close agreement both with each other and with the computed Fay-Riddell two-dimensional 
stagnation point Stanton number. 

times d ThC ^ mPeratUre Pr ° flJe thr ° Ugh thC substrate of one of the gages on the Upilex/Macor cylinder at different 
hmes drnmg the test is shown Figure 20. The discontinuity in the slope of the temperature profile at 0.0508 mm is 

the interface between the low theimal conductivity Upilex and the high (relative to Upilex) conductivity Macor. 


CONCLUDING REMARKS 


thto AF0R r COmPUBr COdC ' IDHEAT ’ deve, °'* d f0r ,be — r dato from 

22:7 27 T7 p,es ' ™ s code ~ S “ “ - — — fto toe 
7 ‘ he conduction problem. The amdydca! techniques are ba*d o„ the conston, 

S “ bS “ e The «»« - -e vafradon of thermal propetries wito 

.emperamre are included through .idler a vafrahle toe™to conductivity, constant thermal diffusivity ppLmadon 

or through Ihe use of an emptoctoly-dedved cofrection factor to the constant .henna, propenie ^ 7 

_ technique is aone-^iensional fruite-volume discredxadon wid, vmialde suhsj dial « " 

fimte-volume method adds the capability for beat d^sfer computes wid, a subsume composed of muldple layern 
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of different materials or of finite wall thickness. Use of this code has been demonstrated with the reduction of 
sample experimental thin-film gage and coaxial thermocouple data sets. 
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APPENDIX A. ANALYTICAL HEAT TRANSFER MODELS 


The governing equation for one-dimensional heat conduction in a semi-infinite solid with constant thermal 
properties can be expressed as (ref. 2) 


d^djxj) 1 cW£ ,0 

dx 1 a 0 dt 


where 


0(x f t)=T(x 9 0-T(x 9 0) 


with boundary conditions 


<?( 0,0 = — * 0 


dx 


= 0 


0(jc,O) = O 


(A-l) 


(A-2) 


(A-3) 


This equation may be solved using Laplace transforms. The transformed differential equation is 


r d 2 e , 

dx 2 


ctoCGrr) = £(^r) 


,de . 

dt 


with boundary conditions 


q(0,s) = -k 0 


d0(0,s) 

dx 


0Ks) = 0 
0 (jc,O) = O 


(A-4) 


(A-5) 


where 


(A-6a) 

£[q(.x,t)] = q(x,s) (A-6b) 
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The Laplace transforms of the derivatives are 


BQ 

£(— ) = s£(0)-0(x,O) 


and 


Now, let 





dx 1 


v=c[e{x,t)] 


then the transformed equation is 


d 2 v 

sv-Q(x,0)= a 0 —tr 
dx 


which has the general solution 


X 

v{x,s) = C[9{x,t)} = Ae 




, fl(x,0) 

s 


Application of the boundary conditions leads to 


/l = 0 


D _ q(0,s)Ja^ 

— n — 

k 0 S 2 

then 





where 


(A-7) 

(A-8) 

(A-9) 
(A- 10) 

(A-ll) 

(A- 12a) 
(A- 12b) 

(A- 13) 
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(A- 14) 


Po = 



- ^Po C p,0^0 


which evaluated at the surface (x=0) gives 

£[<KO,O] = <?(O,5) = /y*£[0(O,O] 

or by the (t=0) boundary condition in (A-3) 


A> 


<7(0, s ) = J jJ { s£ [ 6(0, 0] - 50(0,0)} 


(A- 15) 


(A- 16) 


then 


<7(0,5) = 


Po_ 

V7 


f <0 '° 


(A- 17) 


The Laplace transform of the system can then be inverted through a convolution to obtain an expression for 
surface heat transfer as a function of time 




r' 


» <^(0,A) 

[<7(0,5)] = ?(0,0 = AL [ M ^dX 

n Vt-A 


This equation can be simplified through integration by parts as 


*2 a 2 

| udv = uv )** _ J v du 

Al A] 

where 
« = «(A) 

v= v(A) 


where 


u = 


I 

V <— x 


(A- 18) 


(A- 18) 


(A-20) 


(A-21a) 
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, dX 

du — 

2U-X) 1 

(A-21b) 

, dQ .. 
dv = — — dX 
a IX 

(A-22a) 

* \ 

o 

II 

(A-22b) 


An arbitrary constant can be added to v without changing dv. This is done in order to avoid a singularity later in the 
derivation. 


It 

(A-23a) 

o 

cb 

1 

/— s 

<-< 

o 

Ob 

II 

(A-23b) 


then 


<7(0 = <7(0,0 = 4^ 


0(0,0 , 1 f 0(0,0 - 0(0, A) ^ 
V7 2 J /- 


('-A) 1 


(A-24) 


In order to avoid the singularity at <7=0, Cook and Felderman (ref. 14) reduce the integral to a summation 
by approximating the temperature as a piecewise linear function 


e(0,X) = 9( A) = 0(f M ) + 2M eu ‘-'\ A - ) 

Ar 


(A-25) 


Then by transforming the integral into a summation of discrete integrals 


Vtt . 0 1 


(^-e,_i) f — - (e - g ‘-» } r -Azi- 

, J , On-A) 1 A/ 


f — '\ dXl 

L G--A)* 


(A-26) 


with integration by parts of the second integral 
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(A-27) 


Vtt . 9 n 1 V? I 


1=1 


(dn-Q.-i) 

1 

i 


to 

1 

T 

'W' 

* _f dx 1 


A t 

<i-i 

('„-*) 

IC'n-U* J 


evaluation of the all terms gives 

On-Oj , m-e.-y) 

+ ,)\ 



(A-28) 


The first term in the summation is undefined at (i=n). However, by l'Hopital's rule, the limit as i approaches n of 
this term is zero. 

With the condition that 9(i 0 ) = 0 where t 0 = 0 , the summations can be separated to get 




(A-29) 


Finally, the indices on the second summation are shifted to get 


?o»)= 



-ti+Jtn-ti - 1 


(A-30) 


Equation (A-l) is valid only when the variation in the thermal properties can be neglected, that is, when the 
temperature rise is small. In hypersonic heat transfer experiments the temperature rise is usually quite large, 
therefore, the variation in the thermal properties must be accounted for in the heat transfer model. 

For some materials, such as the metals in a coaxial thermocouple, the thermal conductivity increases 
rapidly with temperature, but the thermal diffusivity increases only slightly. Thus the relationship between heat flux 
and temperature for a coaxial thermocouple can be approximated using the approach suggested by Cook (ref. 3), in 
which the thermal conductivity is allowed to vary, but the diffusivity is held constant. The dependent variable is 
transformed by 



(A-31) 


which leads to the system 
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<? 2 0(x,O _ 1 dQjxj) 
dx 2 a 0 dt 


(A-32) 


The form of equation (A-30) is identical to that of equation (A-l). Thus, it has the solution 

+ y*n _, i-l 


(A-33) 


For insulative materials such as Macor and quartz, which are used as substrates for thin-film gages, 
temperature has a strong effect on both the thermal conductivity and the thermal diffusivity. Thus, neither (A-l) nor 
(A-32) are valid representations of the relationship between heat flux and temperature. It has been found however, 
that accurate correction factors for the effects of variable thermal properties can be determined which can be applied to 
the constant properties solution . The corrected equation takes the form of 


fy(T) =4f3 0 U + P'AT s ] (A-34) 

where the constant property solution can be obtained using the direct method of equation (5) or the indirect method of 
equation (7). 

The correction factors are derived from heat transfer rates computed using the finite-volume technique. The 
finite- volume technique is discussed in Appendix B, while the derivation of the correction factors is presented in 
Appendix C. These corrected semi-infinite substrate heat transfer results are typically within a few percent of the 
results from the more rigorous finite-volume computations. The advantage of this technique is the simplicity with 
which the numerical summations in equations (5,7) can be calculated and then corrected for property variations by 
equation (12). In contrast, the finite-volume solution generally requires more time due to the need to construct a 
finite-volume input file and to analyze the effects of grid spacing on the computations. For this reason it is less 
useful for "real-time" data analysis during wind tunnel testing. However, the finite-volume method is still very 
valuable in that it involves no empiricism, the empirical correction factors for the analytical method are derived from 
it, and it is the only way to correctly solve the finite wall thickness problem or the multiple-layer substrate problem. 
It is thus suggested that the analytical semi-infinite methods be employed to quickly reduce data during a test series, 
and then later that a more rigorous analysis be carried out using the finite-volume technique. 
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APPENDIX B. FINITE- VOLUME HEAT TRANSFER MODEL 


In theory, the most accurate method for computation of heat flux in a solid body would be through the use 
of a three-dimensional finite- volume technique to compute the temperature distribution inside the body in question. 
In this manner, the effects of the variation in thermal properties with temperature, as well as those of lateral and 
transverse heat conduction could be taken into account. However, it would in practice be very difficult to instrument 
a model with enough sensors to accurately measure the entire surface temperature distribution, which would be 
required as a boundary condition for a computational solution. The computational time would also be much greater 
for a three-dimensional solution, and considerable time would be required to generate a three-dimensional 
computational grid and perform a grid resolution study. While accurate, miniaturized temperature sensors and 
powerful computational resources are available today, and thus a three-dimensional computation can be performed, 
this is not a practical approach for a real-time analysis code which can be applied to generic configurations with a 
minimum of preparatory work. For these reasons, the numerical model is presently restricted to that of one- 
dimensional heat conduction. 

An implicit finite- volume discretization is implemented in 1DHEAT. The geometry and notation for this 
model are shown in Figure 21. The temperature distribution within the solid can be obtained from an energy balance 
within each cell: 
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This model is second-order accurate in space and first-order in time. Note that all material properties values 
in the coefficients of T are lagged and evaluated at the (n-7 ) level. Equation (B-3) is valid for the computational cells 
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on the interior of the domain. The temperatures at the exterior and interior (it is assumed that the actual wind tunnel 
model is hollow in order to allow space for the sensor wires) of the domain are functions of the boundary conditions. 

The temperature of the first computational node is extrapolated from the temperature measured on the 
external face of the model by the sensor (Figure 22a.): 

r " = f [ T * age + 2 T ^) (B-4) 


By substitution into equation (B-3), the temperature at node 2 is 
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Several boundary conditions for the back surface are possible: adiabatic surface, constant temperature, and 
measured temperature. The measured temperature boundary condition would be the most accurate, but since it is not 
usually practical to make this measurement, the adiabatic surface boundary condition is generally used. The rationale 
for this boundary condition is that the heat conduction to the stagnant gas in the hollow interior of the model is 
negligible in comparison to the heat flux through the solid. At current, the constant temperature and adiabatic 
surface boundary conditions, but not the measured temperature boundary condition, are implemented in 1DHEAT. 
These boundary conditions are illustrated in Figure 22b. 

The adiabatic boundary condition is implemented by setting to zero in equation (B-l) for the (m-l)th 
computational node 
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For the constant interior surface temperature boundary condition, the temperature at the nuh node is extrapolated 
from the surface by 
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Then by substitution into (B-3) 
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The system of equation (B-3) with boundary conditions of equations (B-5), and (B-6) or (B-8) forms a 
tridiagonal matrix of unknowns which can rapidly be solved using the Thomas algorithm. 

The finite-volume technique in 1DHEAT can be also used to solve the one-dimensional heat conduction 
equation in a substrate with up to three substrate layers in which material thickness, grid spacing, and thermal 
properties are different in each layer. The number of layers can be increased simply by changing the dimensions in 
the code. 

Note that in a multiple-layer problem the thermal conductivities on either side of the interface may be 
discontinuous. Because of this, it is not consistent to extrapolate the cell interface thermal conductivities, ki and k r . 
by averaging the conductivities of the adjacent nodes, that is 


As proof, consider a material in which one layer was a perfect insulator (k = 0) and the other material of the other 
layer had some non-zero conductivity. At the interface of the two materials, the thermal conductivity from equation 
(B-8) would be non-zero, permitting heat transfer into the insulated layer. 

The interface conductivity is properly defined through the series electrical resistor analogy (ref. 6). The heat 
transfer at the interface is given by 
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where R is defined as the conduction resistance 
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In 1DHEAT, the cell wall areas are constant. From the notation in Figure 21, the conduction resistances 
between the nodes and the interface are: 
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(B-12a) 


R m = T 

A*,-, 

2A 




R: = 


1 


Ax 

k — 

2 A, 


Then from equation (B-10) 


Xj-Xj-l 


k =■ 


R i-l + R i 



^,-iAx,_| | k) Ax, 
2A, 2A X 


(B-12b) 


(B-13) 


By defining 
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equation (B-13) becomes 
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With this definition, the interface thermal conductivity in the above example would have the proper value of 
zero. This formulation is also useful in that it permits the node spacing to differ from layer-to-Iayer, although in the 
1DHEAT code it is uniform within a layer. 

The approximate grid Fourier number based on constant material thermal properties is 


<7 = 


PCpAt 

(Ax) 2 


(B-15) 


Because this is an implicit formulation there are no stability restrictions on a, and the grid spacing Ax may be 
adjusted as desired to make trade-offs between solution accuracy and computational time. The grid spacings are varied 
separately by specifying the number of mesh points in each layer in the finite-volume set-up file. It may be 
advisable to adjust the grid spacings so that the Fourier numbers in each layer are identical, although limited 
investigations have not shown any significant effects of varying a from layer to layer. 
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APPENDIX C. MATERIAL THERMAL PROPERTIES 


The accuracy of the material thermal property data is a crucial, but often neglected, factor in the reduction of 
experimental heat transfer data. Since material properties can rarely be considered to remain constant, their 
dependence on temperature must be known, and must be accounted for in the data reduction scheme. The effects of 
thermal properties can be rigorously dealt with by employing a finite-volume data reduction method as detailed in 
Appendix B, or in a simpler empirical manner by the application of correction factors to an analytical constant 
properties solution. However, even if the correct theoretical models are used, it is extremely difficult to obtain 
accurate thermal property data for use in those models. Material handbooks generally only give data at a single point 
(usually room temperature). Thermal property data in the aerospace literature is also lacking, and what data exists is 
often contradictory. 

While useful qualitative information can still be obtained, it is simply not possible to make quantitative 
conclusion from experiments without accurate thermal properties data. This is especially significant in regards to 
CFD code calibration exercises, as discussed by Neumann (ref. 17). Therefore, the thermal properties data 
incorporated into the 1DHEAT code are presented here along with a survey of data sources. It should be noted that 
thermal properties data in 1DHEAT are in some cases different than those published elsewhere. These thermal 
properties represents the best data currently available to the author, but are subject to change and should not be 
considered as THE definitive set of thermal properties data. This data may be modified by the user or data for other 
materials may be added to the 1DHEAT code as discussed at the end of this appendix. Comparisons between 
experimental results based on these properties and computational results are included at the end of this Appendix. 


MATERIAL PROPERTY CURVE FITS 

A summary of room -temperature properties of the materials for which data has been incorporated into 
1DHEAT (chromel, constantan, Macor, quartz, Pyrex, 17-4 stainless steel, and Upilex) is given in Table 4. Curve 
fits for the variations in thermal properties with temperature are given in the respective sections of this appendix for 
each material. These curve fits are all considered to be valid for temperatures up to 600 ”K. 

Type-E Thermocouple 

The Type-E chromel-constantan thermocouple manufactured by Medtherm is composed primarily (80-90% 
by volume) of chromel. Because of this, and because of the uncertainty in constantan thermal properties, it is 
recommended that the chromel curve fits of equation (C-l) be used to represent the Type-E thermocouple. 
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Chrome! 


Because chromel is a metal, it has high thermal conductivity and thus experiences a relatively low 
temperature rise (compared to an insulator) when exposed to heating. The effects of variable material thermal 
properties are thus less significant. However, there is a trade-off in that metals also have higher thermal diffusi vity , 
and thus the time for which the semi-infinite approximation is valid is much shorter. 

Chromel thermal property curve fits are presented in Figures 23a-23e. Chromel specific heat data are taken 
from Touloukian et al (ref. 18), which was the only source found for variation of chromel specific heat with 
temperature. Chromel thermal conductivity data was taken from Touloukian et al (ref. 19), Hoskins 1 (a supplier of 
chromel alloy) and from Medtherm 2 (a manufacturer of the Type-E thermocouples). Thermal diffusivity and thermal 
product values were computed from these data. Material property curve fits based on these data are: 


p = 8714 (kg/m 3 ) (C-la) 

c p = 386.25 + 0.23981 T (J/kg-‘K) (C-lb) 

k = 11.845 + 1.9132 - 10" 2 -T (W/m-’K) (C-lc) 

a = 3.5995 10' 6 + 2. 9656- 10' 9 T - 9.1293- 10 13 T 2 (m 2 /s) (C-ld) 

P = 6398.4 + 6.6331 T (W-s^/m 2 -*K) (C-le) 


Constantan 

Constantan thermal properties were obtained from the same sources as chromel properties. However, less 
data were available in these references and the curve fits below are of higher uncertainty than those for chrome!. 
Constantan thermal properties curve fits are: 


p = 8906 (kg/m 3 ) (C-2a) 

c p = 318.54 + 0. 25745 -T (J/kg-*K) (C-2b) 

k = 8.5591 + 4.6562 -10' 2 T (W/m-’K) (C-2c) 

a = 3. 1797 • 10' 6 + 1. 1966 • 10' 8 T - 4.5508-10' 12 -T 2 (m 2 /s) (C-2d) 

P = 5336.2 + 11.754-T (W-s^/m 2 -‘K) (C-2e) 


These curve fits are plotted in Figures 24a-24e. 

Macor 

Macor is a machinable glass ceramic manufactured by the Corning 3 company. Macor is a thermal 
insulator, which makes it a good substrate material for thin-film gages because insulators have longer semi-infinite 
test times. However, temperature has a strong influence on the thermal properties of Macor, which complicates the 


1 Hoskins Manufacturing Company, Detroit MI 48208 

2 Medtherm Corporation, Huntsville, AL 35804 

3 Corning Incorporated - Advanced Materials Department, Coming, MY 14831. Macor is a trademark of Coming Inc. 
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analysis of heat transfer data. Furthermore, the variation with temperature is not well documented. Thermal 
property data has been assembled from a number of sources, most of which are not in good agreement with each 
other. Macor thermal properties data were obtained from Coming, Miller at the NASA Langley Research Center 
(refs. 5,20), Wannenwetsch et al at the Arnold Engineering Development Center (ref. 21), CALSPAN (cited in ref. 
5), Mentre and Consigny (ref. 22) and Soos (refs. 23-25) and Campell (ref. 26) at the Holometrix 4 company . The 
merits of each data set can be assessed from study of the references cited. Ultimately, Macor thermal property curve 
Fits were based on references 23-25 and 26 because these data are recent (1992-1995), are repeatable, and cover a large 
range of temperatures. 

Macor thermal properties curve fits are : 

p = 2543.84 - 8.0 -lO' 12 (kg/m 3 ) 
c p = 114.04 + 2.5196-T - 1.5136 10' 3 -T 2 (J/kg-’K) 

k = 0.33889 + 7.4682 10' 3 -T - 1.61 18 10 s T 2 + 1.2376- 10' 8 T 3 (W/m-’K) 
a= 1. 3003-10‘ 6 -2.2523-10' 9 T + 1.8571 10 12 T 2 (m 2 /s) 

P = 754.2 + 3.7201 -T - 2.4883- 10' 3 -T 2 (W-s K / m 2 -”K) 

Macor thermal properties are plotted in Figures 25a-25e. 

Quartz 

Like Macor, quartz is a thermal insulator and is thus also a good substrate material for thin-film gages. The 
main drawback to quartz is that it is a more difficult material to machine than Macor. Quartz thermal property data 
were obtained from Coming, GE 5 , NASA (ref. 5) and Holometrix (refs. 23, 26) and are presented in Figures 26a- 
26e. The curve fits presented are based on data from Holometrix for the same reasons as discussed in the section on 
Macor. 

The curve fits for quartz thermal properties are: 

p = 2192.5 (kg/m 3 ) 

c p = 197.85 + 1.9893 - T - 7.3896-10" 4 -T 2 (J/kg-’K) 
k = 0.96157 + 9.5491 10 -4 -T + 5. 5465 - 10' 7 - T 2 (W/m-°K) 
a = 1.5191 10* 6 -4.136 10* 9 T + 7.2707- 10' !2 -T 2 - 4.4242 -10' 15 -T 3 (m 2 /s) 

P = 805.48 + 2.1192-T (W-s M /m 2 -°K) 


(C4a) 

(C-4b) 

(C-4c) 

(C4d) 

(C4e) 


(C-3a) 

(C-3b) 

(C-3c) 

(C-3d) 

(C-3e) 


4 Holometrix Incorporated, Bedford, MA 01730-2323 

5 GE Quartz Products, Cleveland, OH 44117 
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Pyrex 


Like quartz and Macor, Pyrex 6 is another insulative material used as a thin-film gage substrate. The curve 
fits presented here are from reference 5 and are presented in Figures 27a-27e. Because Pyrex thermal properties vary 

much more with temperature than those of other materials, and it is not easily machined, it is generally not a 
preferred substrate choice. 


p = 2227.0 (kg/m 3 ) 

Cp = -85.414 + 4.2839 T - 5.7482 - 10‘ 3 -T 2 + 3.1047 -Hr 6 T 3 (J/kg-’K) 
k = 1.5146 - 5.9068 ■ 10' 3 -T + 1.8165 • 10 s -T 2 (W/m-'K) 

« = 1.0622- 10' 6 -3.0254 10‘ 9 ■ T + 7.3053 • 10 12 T 2 (m 2 /s) 

P = 541.63 + 2.0121-T + 4.3415-10' 3 T 2 (W-s^ /m 2 -°K) 


(C-5a) 

(C-5b) 

(C-5c) 

(C-5d) 

(C-5e) 


Upilex 

A technique currently being developed at NASA Langley is the application of thin-film gages to a layer of 
Upilex film which is then applied to a wind tunnel model. Upilex is a thin (1 to 5 mil thickness) polyamide film 
manufactured by Ube 7 which has very low thermal conductivity and a highly finished surface. Because of the surface 
quality of the film, it is possible to obtain higher gage quality and spatial resolution than with other substrate 
materials. However because the Upilex is so thin, it usually cannot be treated as a semi-infinite substrate, and thus 
the finite-volume method must be used to carry out a multiple-layer analysis which takes into account both the 
Upilex surface layer and the sub-layer (typically Macor). Because this is a relatively new material, there is little data 
available on the thermal properties. Curve fits given below are based on data from Ube and Holometrix (ref 26) and 
should be considered as preliminary. These curve fits are plotted in Figures 28a-28e 


p = 1490.0 (kg/m 3 ) 

c p = -2258.2 + 19.492 -T - 3.7267 -1 O' 2 -T 2 + 2.7812 -10‘ 5 T 3 (J/kg-‘K) 
k = -0.26918 + 3.6348 10' 3 T - 7.2432 10' 6 -T 2 + 5.0056 1 O’ 9 T 3 (W/m-TC) 
a = 5.9629 10' 7 -2.2994 10' 9 T + 4.2071 10' 12 -T 2 - 2.801 1 - 10 15 T 2 (m 2 /s) 

P = -1068.7 + 10.619- T - 2.0458 • 10' 2 -T 2 + 1.4544 10 s -T 3 (W-s^ /m 2 -TC) 

17-4 Stainless Steel 

A technique sometimes utilized for wind tunnel model construction is to fabricate a model from steel, 
aluminum, or other easily machinable material. Sections are then cut out from the model and replaced with 
instrumented inserts. These inserts are typically Macor, quartz or Pyrex substrates with thin-film gages. Ideally, the 
substrates are thick enough that they can be considered as semi-infinite. However, this problem can be analyzed as a 


6 Pyrex is a trademark of Corning Inc. 

7 Ube Industries, New York, NY 10103 


(C-6a) 

(C-5b) 

(C-6c) 

(C-6d) 

(C-6e) 
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multiple-layer problem using the finite-volume technique. Theimal property data for 17-4 stainless steel obtained 
from ARMCO 8 is thus given here for use in a multiple-layer analysis. However, this data should only be taken as 

an example of metallic substrate properties since there are many varieties of steel or other metals which can be 
employed in model construction. 


p = 7800 (kg/m 3 ) 
c p = 460 (J/kg-'K) 

k = 11.63 + 1.4816 10' 2 T (W/m-*K) 
a = 3.2414 10' 6 +4.1293 10’ 9 -T (m 2 /s) 
= 6601.2 + 3.3208 T (W-s K /m 2 -”K) 


(C-7a) 

(C-7b) 

(C-7c) 

(C-7d) 

(C-7e) 


CORRECTION FACTORS 


The numerical method discussed in Appendix B can be used to obtain an "exact" solution to the problem of 
one-dimensional transient heat conduction in a solid with variable thermal properties. However, this is not 
considered to be the optimal method for real-time reduction of large amounts of data because the finite-volume 
method generally requires more set-up work and analysis than the analytical method. This makes it somewhat less 
user-friendly, because non-experimental factors such as grid resolution and boundary conditions must be considered 
before the 1DHEAT code can be used to reduce data. The preferred methods for real-time data reduction are the semi- 
infinite substrate methods of equations (5) or (7) with the correction for variation of material properties with 
temperature of equation (12) for thin-film gages, and the semi-infinite substrate methods of equations (9) or (10) for 

coaxial thermocouples. The finite-volume method should be used to reduce selected sets of data to verify these 
results. 

Correction factors for Macor were derived in a manner similar to Cook (ref. 3) and have the form 

W) = fyo f 1 + P' AT s ] (C-8) 

The correction factor, 0 , for Macor was determined by comparing the computed heat transfer rates for a parabolic 
temperature-time history, T = T(^), (corresponding to a constant heat transfer rate for a substrate with constant 
thermal properties) for a substrate with thermal properties held constant to a substrate with variable thermal 
properties. In both cases the heat transfer rates were computed using the finite-volume method. These heat transfer 
rates are plotted versus temperature rise in Figure 29. The correction factor shown can be derived from this figure as 


8 Armco Stainless Steel Products, Baltimore, MD 21203 
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(C-9) 


PUT,) 



AT-, 


The correction factor is considered to be valid for temperatures up to 600 *K and is implemented in 1DHEAT using 
the linear fit shown. 


P' = 7. 380 • KT 4 - 4. 604 ■ 10" 7 • AT, 


(C-10) 


Note that although the linear fit is less accurate at lower temperatures (the computed correction factor goes to zero at 
room temperature) the actual correction to the heat transfer rate is the correction factor times the temperature rise, and 
the error is thus negligible at low temperatures. 

The finite-volume solutions for constant and variable thermal properties and the corrected analytical solution 
are plotted in Figure 30 as a function of time. For this parabolic time history, the agreement is excellent. Although 
the correction factor values are derived only from this parabolic temperature-time history, they appear to be valid for 
any temperature-time relationship within the temperature range given for the material property curve fits. A sample 
comparison between corrected analytical results and finite-volume results is shown in Figure 31 (stagnation point 
heating from the demothin" data set) and the agreement is excellent. 

The correction factor for quartz was derived in the same manner as for Macor. The thermal properties for 
Pyrex in (C-5) are taken from reference 5, so the same correction factor is used. 

For quartz, the correction factor is 


P = 9.414 • 10 -4 - 8.018 - 10 -8 AT, (C . j 1} 

The correction factor for Pyrex is a constant: 

P = 2.33 -KT 3 (C-12) 

No correction factors were determined for the metallic materials or for Upilex. For the metals, the surface 
temperature rises experienced are generally small enough that the constant diffusivity approximation (eq. 7) is valid, 
and this method is implemented in the 1DHEAT code. Data from thin-film gages on Upilex are analyzed by the 
finite-volume method so no correction factors are required. 
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ADDITIONAL MATERIALS 


Users wishing to incorporate new materials into the 1DHEAT code or modify current data must supply 
curve fit coefficients for the thermal conductivity, theimal diffusely, thermal product and thermal product correction 

factors. These additions must be made in the "matprops" subroutine, and the choices must be added to the menu in 
the "readsetup" subroutine 

The thermal diffusivity and thermal conductivity values for a material are defined in 1DHEAT by fourth- 
order polynomial curve fits (for lower-order curve fits, elements of the coefficient arrays may be set zero). 


Thermal conductivity: 

k(T) = a[ k) +4 k) -T + a W ■ T 2 + a {k) ■ T 3 


Thermal diffusivity: 

a(T) = a[ a) + • T + - T 2 + a[ a) ■ T 3 


(C-13) 


(C-14) 


The thermal product and correction factor are defined by a three element array, /?, , where /J, , is the room 

temperature value of the thermal product, and the two remaining elements are the coefficients of a curve fit for the 
correction factor 


P -Pl+PlW-Tamb) (C-5; 

For metallic substrates that are represented by the variable thermal conductivity model, the first element is 

again the room temperature value of the thermal product, and the second and third elements are the coefficients for a 
thermal conductivity curve fit: 


T~ = Pz + Pi'T 

K o 


(C-6) 


VALIDATION OF THERMAL PROPERTIES DATA 

As mentioned previously, the thermal properties data given in this work are in some cases quite different 
from data presented elsewhere. These differences are most pronounced for the materials Macor and quartz. Use of the 
current properties is validated by the results presented in Figure 32. This figure show the Stanton number 
distributions on two 1" radius hemisphere models tested side-by-side in the NASA LaRC 31" Mach 10 Air Tunnel. 
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One hemisphere was machined from quartz and one from Macor, and both were instrumented with thin-film 
resistance gages. These distributions are compared to a numerical solution computed using the viscous shock layer 
technique (courtesy of Roop Gupta at NASA Langley) and all distributions are normalized by the computed Fay- 
Riddell (ref. 27) stagnation point Stanton number. Both the quartz and Macor results obtained using the current 
values are within 5% of the computed VSL solution. It was found that experimental results computed using thermal 
property data from other references varied between 10% to 20% from the VSL solution. 
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APPENDIX D. 1DHEAT INPUT FILES 


SETUP FILE "lDdemocoax.inp" 


Listing of sample setup file " lDdemocoax.inp" 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


Sample Coaxial Thermocouple Hemisphere Data 


300 .2d0 
1419 . 7d0 
4 . 3 078d-3 
0 . 7603d6 
300 . OOdO 
Air 
2 . 5d0 
0 . 02d0 
300 
4 . OdO 
4 . 5d0 
l.OdO 
2 

S/R 


ambient temperature (deg K) 
freestream velocity (m/s) 
freestream density (kg/m 3 ) 
total enthalpy (h_tot - h_298) (J/kg) 
reference temperature (deg K) 
test gas (Air, He, N2 , CF4 , C02 ) 
time zero (sec) 
data sampling period (sec) 
number of points in integration window 
start time of averaging window (sec) 
end time of averaging window (sec) 
adiabatic wall recovery factor 
columns of position data 
title of position column #1 


16 

theta 

title 

of position 

column #2 


17 

GAGE ID 

S/R 

theta 

ALPHA_R 

GAGE TYPE 

18 

' TC7 1 

-1.5709 

-90.0 

0.00 

coax 

19 

1 TC6 * 

-1.3091 

-75.0 

0.00 

coax 

20 

1 TC5 ' 

-1.0472 

-60.0 

0.00 

coax 

21 

1 TC4 ' 

-0.7854 

-45.0 

0.00 

coax 

22 

1 TC3 ’ 

-0.5236 

-30.0 

0.00 

coax 

23 

' TC2 1 

-0.2618 

-15.0 

0.00 

coax 

24 

’TCI • 

-0.0000 

0.0 

0.00 

coax 


25 


end 


Description 

of "lDdemocoax.inp" file 

Line 1: 


Text header describing the data. 

Line 2: 

T 

amb 

Ambient temperature before the run 

Line 3: 


freestream velocity 

Line 4: 


freestream density 

Line 5: 

A/W 

total enthalpy, A h wl = A(r w )- *(7^) 

Line 6: 

t,/ 

reference temperature at which to extrapolate heating rates 

Line 7: 

gas 

identifies test gas. Must be "Air", He”, "N2”, "CF4" or M C02' 

Line 8: 

*0 

time zero. Time at which data reduction begins 

Line 9: 

At 

time between data samples. 

Line 10: 

n ,o, 

number of data points after time zero to continue data reduction 

Line 11: 

h 

start time of window over which reduced data is to be averaged 

Line 12: 

h 

end time of window over which reduced data is to be averaged 


SUBSTRATE 
chrome 1 
chrome 1 
chrome 1 
chrome 1 
chrome 1 
chrome 1 
chrome 1 
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Line 13: 

Line 14: 

Lines 15 -16: 
Line 17: 

Line 18-24: 


r adiabatic wall recovery factor used in Stanton number calculation 
n col number of columns (maximum of 4) of position data to be read 

titles of position data columns. 

Text header identifying data in following lines. 

Gage data lines which contain: 

gageid(n): identifier for nth gage, may be up to 6 characters 

pos(n,k): numerical position data for nth gage corresponding to kth position column 

data title 

a R ( n ) temperature coefficient of resistance (in l/'R) for the nth thin-film gage. 

This is a dummy value for coaxial thermocouples or if the raw data is 
already converted to 'K 

gage type either "coax" (for coaxial thermocouple) or "thinfilm" (for thin-film gage) 

substrate substrate material. For multi-layer substrates, this identifies the material 


properties to be used to reduce the data if the analytical method is chosen and the 
substrate is approximated as a single layer 
Line 25 Terminates reading of gage data 


The flow properties information in the setup file is not required for the computation of the heat transfer 
rates. However, it is used in the computation of the Stanton number, C h , the heat transfer coefficient, H, and the 
reference heating value, q re j. It must be emphasized that the heat transfer coefficient computations are only as 
accurate as the user-input flow properties information. Furthermore, one of the inherent assumptions in these 
computations is that the flow is steady; that is, the density, velocity and total enthalpy do not vary with time. For 
unsteady flow, these results will not be valid. 

Wall enthalpy as a function of temperature is also required for the heat transfer coefficient computations. 
Currently, 1DHEAT incorporates enthalpy curve fits for air, N 2 , He, CF 4 and C0 2 gases. The curve fit coefficients 
are taken from McBride (ref. 28), and are valid for wall temperatures between 300°K and 1000‘K. If needed, 
additional curve fits may be added in the "hwall" function of the code, and the additional choices can be added to the 
menus in the "readsetup" subroutine. 

Several times are required in the input file. The times t 2 and t 2 define the start and end of the window over 
which the data is to be averaged. Initial selection of the proper averaging window in the input file is not critical, as 
the data can easily be re-averaged using 1DHEAT code options 3 or 7. In fact, a better definition of the averaging 
window can be made by first running the code with assumed window start and end times, and then analyzing the heat 
transfer-time histories to determine the optimum values for tj and t 2 . 

The time t 0 defines the point at which reduction of the data begins. The exact definition of t 0 is also not 
critical as long as t 0 is before the gages experiences any heating. In fact, there is a built-in safety factor on t 0 in this 
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code that will step back 50 time increments from the user-input value (or fewer if there are less than 50 data points) 
Note that the 1DHEAT code will compute baseline pre-test voltages (for thin-film gages) for each sensor by 
averaging over the first 25 points from t 0 . The 50-point safety factor and 25-point baseline averaging values are set 
in the 'setup' subroutine of the code. 

Entries in the input file must correspond to allowable values. Allowable text strings for the test gas are 
"Air", "N2", "CF4" "He" and "C02". The gage type can be "thinfilm" for thin-film resistance gages or "coax” for 
coaxial surface thermocouples. Finally, the substrate type can be set to "macor”, "quartz", "pyrex", "chromd", 
"constantan", "stainless", or "upilex". 

FINITE-VOLUME SETUP FILE "lDdemocoax.fvinp" 

Listing of sample setup file "lDdemocoax.fvinp" 

1 SAMPLE COAX DATA FINITE -VOLUME SETUP FILE 


2 

GAGE ID 

LAYERS 

WALL BC 


3 

Layer # 

NODES 

THICKNESS 

SUBSTRATE 

4 

1 TC7 ' 

1 

0 


5 

1 

43 

1 . 216.-2 

chrome 1 

6 

' TC6 1 

1 

0 


7 

1 

43 

1 . 27d-2 

chrome 1 

8 

1 TC5 1 

1 

0 


9 

1 

43 

1 . 27d-2 

chrome 1 

10 

' TC4 ' 

1 

0 


11 

1 

43 

1 .27d-2 

chrome 1 

12 

1 TC3 1 

1 

0 


13 

1 

43 

1 . 27d-2 

chrome 1 

14 

' TC2 ' 

1 

0 


15 

1 

43 

1 . 27d-2 

chrome 1 

16 

'TCI ’ 

1 

0 


17 

18 

1 

end 

43 

1 . 27d-2 

chrome 1 


Description of "lDdemocoax.fvinp" file 

Line 1: Text header describing the file 
Line 2: Column titles 
Line 3: Column titles continued 
Lines (4,5), (6,7), etc.: 

The first line in each pair contains 

gageid(n): identifier for nth gage, may be up to 6 characters 

nlay(n): number of substrate layers of nth gage 

nbc(n): interior boundary condition for nth gage 

(0 for adiabatic, 1 for constant ambient temperature) 
The second line in each pair contains 
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k- layer number 

nplay(njc): number of points in kth layer of nth gage 

d(njc): thickness of kth layer of nth gage 

subs(njc): substrate material of kth layer of nth gage 


Note that in this example, the substrate has only one layer. A substrate with more layers would have 
additional lines of the same form as line 5. For example a thin-film gage on a 2 mil (0.0508 mm) layer of Upilex 

over a 1/4" (0.635 cm) thick Macor layer with 5 nodes in the upilex layer and 244 nodes in the Macor layer would 
have lines of the following form for each gage: 


1 TF1 1 2 

1 -x 5 

2 244 


0 

5 . 08d-5 
6 . 35d-3 


upilex 

macor 


INPUT TEMPERATURE FILE " lDdemocoax.degk" 


Listing of sample input I 

1 Time TC7 TC6 

2 0.0000 300.25 300 

3 0.0200 300.24 300 

4 0.0400 300.21 300 

5 ... (remainder of fill 


" lDdemocoax.degk" 



TC5 

TC4 

19 

300.23 

300.21 

19 

300.24 

300.19 

18 

300.21 

300.18 

omitted) 



TC3 TC2 TCI 

300.21 300.23 300.22 

300.21 300.23 300.22 

300.19 300.21 300.21 


Description of "lDdemocoax.degk" file 

Line 1: Text header describing following lines. 

Lines 2+: Contain the time index followed by the temperatures (or voltages in a ".volt" file) 

for each gage at that time index. Reading is terminated by the end of the file 


Note that a voltage input file (with a '.volt' extension) will have the same format except that temperatures 
are replaced with voltages. 
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TABLE1. 1DHEAT data files 
File template: lD testlDrunlD .exr 


Extension Description 


.inp 

- setup file 

.fvinp 

- finite- volume setup data 

: volt 

- raw data (volts) 

.degk 

- raw data (°K) 

A 

- time history of sensor temperatures 

qx 

- time history of heat transfer 

.chx 

- time history of Stanton number 

.qresx 

- time history of heat transfer residual 

.chresx 

time history of Stanton number residual 

.distx 

- time-averaged heat transfer data 

. distxnew 

- recomputed time-averaged heat transfer data 

Aback 

- back face temperature-time histories 

.gageid 

- temperature profile though substrate of sensor 
"gageid" 


x is the data reduction method: 

1) indirect semi-infinite 

2) direct semi -infinite 

3) finite-volume 


TABLE 2. Thermocouple Coefficients 
from National Bureau of Standards 


TABLE 3. 31” Mach 10 Air Flow Conditions 


TABLE 4. Approximate room temperature material properties 


MATERIAL 

P 

(kg/m 3 ) 

B88B 

k 

(WAn-*K) 

a 

(m 2 /s) 

P 

(W-s I/2 /m 2 -*K) 

j CHROMEL 

8714 

458 

17.5 

4.40 xlO* 6 

8376 

CONSTANTAN 

8906 

395 

22.4 


8840 

| MACOR 

2568 

731 

1.46 

7.94 x 10- 7 

1642 

QUARTZ 

2192 

726 

1.30 

8.15 x 10- 7 

1435 

PYREX 

2227 

458 

1.37 

■'frrnrgM 

1527 

UPILEX 

1490 

967 

0.303 

2.11 x 10- 7 

661 

17-4 STAINLESS 

7800 

460 

16.0 

4.47 x IQ 6 

7587 



"demothin" 

"democoax” 

Pi (Pa) 

T, (-K) 
u, (m/s) 

Mj 

Re t (1/m) 
Po .2 (Pa) 
T 0 ^ CK) 
Aho. 2 (I/kg) 

65.18 

52.73 

1419 

9.75 

1.570 x 10 6 
8048.8 
1012.0 
0.760 xlO 6 

70.11 

53.04 

1413 

9.68 

1.661 x 10 6 

8528.0 
1004.2 
0.751 x 10 6 


Type-E Thermocouple 


Jk) 


0.104967248 


17189.45282 


a 2 


-282639.0850 


a 3 


12695339.5 


24 


-448703084.6 


a 5 


1.10886 x 10 10 


26 


-1.76807 x 10 11 


a 7 


1.71842 x 10 12 


2 S 


-9.19278 xlO 12 


2q 


2.06132 x 10 13 









































































"1DHEAT" v2 . 20 PROGRAM OPTIONS: 

SEMI -INFINITE METHOD 

1. Reduce raw data (volts). 

2. Reduce raw data (degrees K) . 

3 . Recompute statistics over new 
time interval. 

4. Compute residuals. 

FINITE-VOLUME METHOD 

5. Reduce raw data (volts). 

6. Reduce raw data (degrees K) . 

7. Recompute statistics over new 
time interval. 

8. Compute residuals. 

9. Temperature profile for a gage at a 
specified time. 

10. Quit. 


FIGURE 4: 1DHEAT main menu 
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Figure 24: (a-e) Constantan thermal properties 







MACOR SPECIFIC HEAT 


MACOR THERMAL CONDUCTIVITY 



Figure 25: (a-e) Macor thermal properties 













PYREX SPECIFIC HEAT 


PYREX THERMAL CONDUCTIVITY 



Figure 27: (a-e) Pyrex thermal properties 
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